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Abstract 

Uganda's forests are globally important for their conservation values but are under pressure from increasing human 
population and consumption. In this study, we examine how conversion of natural forest affects soil bacterial and fungal 
communities. Comparisons in paired natural forest and human-converted sites among four locations indicated that natural 
forest soils consistently had higher pH, organic carbon, nitrogen, and calcium, although variation among sites was large. 
Despite these differences, no effect on the diversity of dominant taxa for either bacterial or fungal communities was 
detected, using polymerase chain reaction-denaturing gradient gel electrophoresis (PCR-DGGE). Composition of fungal 
communities did generally appear different in converted sites, but surprisingly, we did not observe a consistent pattern 
among sites. The spatial distribution of some taxa and community composition was associated with soil pH, organic carbon, 
phosphorus and sodium, suggesting that changes in soil communities were nuanced and require more robust 
metagenomic methods to understand the various components of the community. Given the close geographic proximity of 
the paired sampling sites, the similarity between natural and converted sites might be due to continued dispersal between 
treatments. Fungal communities showed greater environmental differentiation than bacterial communities, particularly 
according to soil pH. We detected biotic homogenization in converted ecosystems and substantial contribution of (3- 
diversity to total diversity, indicating considerable geographic structure in soil biota in these forest communities. Overall, 
our results suggest that soil microbial communities are relatively resilient to forest conversion and despite a substantial and 
consistent change in the soil environment, the effects of conversion differed widely among sites. The substantial difference 
in soil chemistry, with generally lower nutrient quantity in converted sites, does bring into question, how long this resilience 
will last. 
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Introduction 

Tropical rainforests (TRF) possess most of the world's terrestrial 
biodiversity and deforestation is the leading cause of biodiversity 
loss [1,2]. Due to their high biodiversity and endemism, the 
tropical rainforests in Uganda's Nile river watershed are among 
the world's most important for their conservation values. But these 
areas are under pressure. The United Nations Population Division 
[3] predicts that the population of the Nile Basin states will 
increase by 57% from 2010 to 2030, reaching 647 million people. 
This rapid population growth, high levels of poverty and prevalent 
civil insecurity continue to exert severe pressure on natural 
resources in the region. Uganda in particular has one of the 
world's highest population growth rates (3.2% per year) [4]. Most 
of this growing population (nearly 80%) is dependent on 



agriculture leading to large scale and continuing conversion of 
natural habitats [5]. 

Soil communities form the foundation of any ecosystem, in 
terms of nutrient cycling and availability, so understanding how 
land conversion affects these communities is an important first 
step. The effect of land use change on soil microbial communities 
has been studied in South American and Southeast Asian forests 
[6,7], but not in the biodiversity hotspots of the Nile river 
watershed. There is considerable global concern about the loss of 
biodiversity and the consequences for human well-being [8]. 
Microorganisms in particular play a vital role in many ecological 
processes and environmental services [9]: these roles are not 
always apparent or well characterized but if all microbes died the 
world would rapidly become buried in undecomposed dead 
material. Due to their significance in maintaining ecosystem 
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function and productivity [9,10], our study ofTers a vital 
exploratory appraisal of microbial community dynamics in natural 
TRF and human-converted sites. We don't know if there are 

reasons to be concerned unless we look. Developing such 
knowledge is critical at this point, because populations in the Nile 
river watershed are highly dependent on forests for basic 
requirements such as food and fuel wood, with the environment 
contributing between 40-60% of the gross domestic product 
(GDP) of the Nile riparian states [11]. 

Because of widespread loss of biodiversit)', focus from species 
conservation within particular habitats has been shifted to 
conservation of communities [12,13]. It is therefore important to 
explore and understand how composition and diversity changes 
across spatial scales in a given context [14-16]. Changes in 
ecosystems caused by conversion to intensive management can 
lead to biotic homogenization, the increase in community 
similarity over time and/ or space and an imphed loss of rare 
and vulnerable taxa when examined at larger scales [17-19]. 
Because microorganisms are the most diverse organisms on earth 
with most taxa and respective functions and behaviors as y(-t 
unknown, determining their sensitivities and biogeography 
remains a major challenge. But in the longer term such knowledge 
will help us better understand the sustainabihty of land-use systems 
and associated environmental vedues. 

This study was therefore necessary as a first step in exploring 
these relationships, and to enhance understanding so as to 
contribute to the informed and appropriate stewardship of the 
region's natural resources. Our objective was to estabhsh how 
forest conversion and soil factors affect soil bacterial and fungal 
diversity and community composition in the tropical rain-forests in 
the Nile river watershed of Uganda. We chose four forest sites 
found within protected areas, with paired treatments within each 
forest; (1) natural and (2) converted ecosystem sites. The natural 
forest ecosystem at each site had suffered minimal human 
disturbance, while converted areas had been transformed to 
cropland. These matched sites found in different locations and 
environmental conditions each experienced different land use 
histories, conservation circumstances and individual challenges for 
management. 

In each matched set of natural and converted sites, we 
compared soil physical and chemical properties and microbial 
community diversity and composition using standard PCR-based 
genotyping techniques. We then calculated community similarity 
indices between sites. This approach would allow us to examine 
both environmental and biotic changes in the soil community 
associat(^d with conversion. Disturbances of sufficient magnitude 
or duration may alter an ecosystem and force a difierent regime of 
predominant processes and structures that favor some populations 
over others [20]. 

We tested the nuU hypotheses that there was no diflference in soil 
properties, band-types, and diversity between treatments. The 
influence of soil properties on microbial community diversity was 
measured by discriminant analysis and canonical correspondence 
analysis (CCA), [21,22]. Because additive partitioning of diversity 
provides a useful framework for quantifying the spatial patterns of 
diversity across hierarchical spatial scales [23,24], we partitioned 
total diversity (y) in each ecosystem type (natural and converted) 
into additive components representing within-community diversity 
(a) and between-community diversity (fi). Our objective was to 
identify the most important sources of total diversity so as to 
propose conservation measures for microbial communities in the 
TRF ecosystems of the Nile river watershed of Uganda. 



Methods 

Site description 

We selected four tropical rainforest (TRF) sites because of their 
relative size, biodiversity, socio-economic and scientific impor- 
tance (Fig 1). Mabira forest is located between the highly 
populated and urbanized Kampala city on the western side; the 
extensive and mechanized Lugazi sugar and tea plantations on the 
Eastern; and Lake Victoria on the southern side. Budongo forest is 
located next to the extensive Kinyala sugar plantations on one side 
and a densely populated mainly subsistence population scattered 
around it. Maramagambo and Kaniyo Pabidi are located within 
Queen Elizabeth and Murchison Falls national parks (NP) 
respectively. These two NP forests had perhaps the best protection 
due to presence of Uganda Wildlife Authority (UWA) personnel. 
However, Maramagambo's location starting on the steep slopes of 
the rift valley subjected it to frequent storms with strong runoff 
flow that swept away most of its top soil (Table 1). 

Soil sampling design 

We collected 400 core soil samples within 40 plots (1000 m^ 
each) in four TRF sites (Fig 1). We sampled five plots from each 
site of the natural TRF and five plots from the converted TRF. We 
established the plots at least 100 m from the ecosystem edge and 
500 m apart and collected 10 evenly placed core subsamples of top 
soil (0-15 cm) from each plot and homogenized them into one 
sample per plot. We then derived a 500 g composite sample from 
the mixture, sieved and packed it for physical and chemical 
analyses and DNA extraction. 

Sample preparation 

We sieved 100 g of the soil on-site through a 4 mm mesh, 
transported it to the laboratory on ice, and stored in a freezer at — 
40°C prior to nucleic acid extraction and analysis. We kept the rest 
of the soil for drying and physical and chemical analysis. We 
performed DNA extractions from 1 g of soil using the Ultra Clean 
soil DNA kit (Mo Bio Labs, Solana Beach, CA, USA) following the 
manufacturer's protocol. The purified DNA was detected by 
agarose gel electrophoresis, and the DNA was amplified by 
polymerase chain reaction (PGR). 

Soil property analyses 

We measured the soil [)11 in 2.5:1 water to soil suspension using 
a pH meter (10 g soil+ 25 ml of distilled water, shaken for 30 min 
and read on a calibrated pH meter). We then used the Walkley 
and Black method [25] to analyze soil organic carbon (SOC) and 
the Kjeldahl method [26] to determine soil nitrogen. We 
measured the soil phosphorus by tlu; Bray and Kurtz no. 1 
method [27]. The photoelectric flame photometer was used to 
determine the soil potassium, sodium and calcium after extraction 
with neutral ammonium acetate. We used the atomic absorption 
spectrometer to measure the soil magnesium after extraction with 
neutral ammonium acetate. The Bouyoucos hydrometer method 
adopted from Gee and Bauder [28] was used to determine soil 
texture. The soil copper and iron were then determined using the 
atomic absorption spectrometer after extraction with EDTA. 

PCR amplification and DGGE analysis 

Polymerase chain reaction-denaturing gradient gel electropho- 
resis (PCR-DGGE) method has been used extensively in microbial 
ecology and is a robust and cost effective method for exploratory 
classification of microbial communities [29]. Following soil DNA 
extraction, we performed a PGR for each DNA extraction to 
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amplify the 16S rRNA genes for bacteria and 18S rRNA genes for 
fungi using universal primers (Table 2). 

PGR reactions had a final volume of 25 |il containing a final 
concentration of 1 x TaKaRa ExTaq PGR buffer with MgGl2, 
300 pM of primers for bacteria. We then added 200 |a,M dNTPs, 
2.5 U ExTaq DNA polymerase (TaKaRa Bio, Otsu, Japan) and 
mUliQ, H2O to complete the volume, BSA was also added for the 
fungal community analysis. We performed PGR cycles with an 
initial denaturing temperature of 95 °G for 5 min, followed by 35 
cycles of 95 °G for 30 sec, annealing temperature of 50 °G for 
30 sec, extension of 72 °G for 1 min; and a final extension of 72 
°G for 10 min. We checked the product of the PGR-rounds and 
quantified by agarose gel-electrophoresis. 

We then performed 16S rRNA and 18S rRNA-DGGE analysis 
using a universal mutation detection system (Dcode Bio-Rad, 
Richmond, GA, USA) with a 6% and 8% acrylamide gel for 
bacteria and fungi respectively containing a gradient of 40-60% 
denaturant (100% denaturant contains 7 mol urea and 40% 
formamide). We applied 100 ng of PGR samples to the DGGE 
gel. DGGE was performed in 1 x TAE Buffer (40 mol Tris/ 
acetate, pH 8; 1 mol ethylene diaminetetra acetic acid) at 60 °G 
and a constant voltage of 60 V for 1 6 hours. After staining with 
SYBR Green 1, we recorded the DGGE gels as digital images and 
analyzed the DNA band numbers using image-processing software 
after subtracting background noise. 

Data analysis 

We used the Rolling disk method with Qjaantity One (Bio-Rad 
laboratories Inc.), which normalizes the band pattern from 
electrophoresis for identification of each band. We then converted 
the band patterns into binary data based on the presence or 
absence of each band for part of our analysis. The DGGE 
fingerprints were interpreted in terms of band richness (number of 
predominant DGGE bands/population). The pixel intensity of 
each band was detected by Quantity One software and is 



expressed as relative abundance (P,) [30]. Shannon Index [H') 
and Simpson index (Z)), the most widely used diversity indices were 
then calculated using the richness and relative abundance data 
following the equations: 

R 

H'=-Y,PiInP. (1) 

!=1 



R 

D=Y.^Pif (2) 

1=1 

Where R, the richness, is number of different bands each data 

set contains, P, = and «, is the abundance of the ilh band and 

A' the total abundance of all bands in the sample. 

Band-type data of the DGGE fmgerprints was then used to 
derive the alpha diversity (bands per sample and ecosystem type), 
beta diversity (total bands per site) [31]. Jaccard's similarity indices 
[32] between converted and natural TRF sites were determined 
using the equation: 

Jaccard's Similarity Index = ^/(^ + _S + C) 
Where, 

A = Total number of bands present in both converted (G) and 
natural (N) ecosystem samples (plots) (also fi-diversity) 
B = Number of bands present in G but not in N 
C = Number of bands present in N but not in G 
We determined the influence of site factors as revealed by soil 
physicochemical properties on the variation of soil microbial 
communities by applying discriminant analysis using Statistical 
Package for the Social Sciences (SPSS). This was done to assess the 
relative importance of each predictor variable (pH, SOG, N, P, K, 
Na, Ca, Mg, and soU texture). We also used the Mann-Whitney 
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Table 2. Sequences of primers used in study. 





Microorganism 


Primer 


Sequence (5'-3') 


Reference 


Bacteria 


F357 


CGC CCG CCG CGC GCG GCG GGC GGG GCG GGG GCA CGG GGG GCC TAG GGG AGG GAG GAG 


[72] 




907R 


CCG TCA ATT GMT TTG AGT TT 




Fungi 


FF390 


CGA TAA CGA AGG AGA GGT 


[73] 




FR1GC 


AlC GAT TCA ATC GGT AIT 





doi:l 0.1 371 /journal.pone.Ol 0481 8.t002 



test to examine difTerences between soil properties in natural and 
converted ecosystems, and microbial communities in natural and 
converted ecosystems. 

We tested the nuU hypothesis that diversity is uniform at all 
spatial scales by additive partitioning of total diversity (y diversity). 
To determine contributions of a and P diversity to overall diversity 
across a range of spatial scales [14,23,33], an additive relationship 
between diversity components (i.e., P = y - ot) was derived (Fig 2.). 
The scale at which diversity is maximized was therefore identified 
[23,34], to facilitate planning processes and management strate- 
gies to conserve natural levels of diversity accordingly [35-38] . 

We used PARTITION 3.0 software [39] to calculate average 
diversity at each scale and diversity was measured as band 
richness. Individual-based randomization procedure in the soft- 
ware was used to test whether the observed partitions of diversity 
within the ecosystem could have been obtained by a random 
allocation of lower-level samples nested among higher-level 
samples [34]. NuU values of |3; obtained from 1,000 randomiza- 
tions were used to obtain a jf)-value for the observed P; at each 
hierarchical scale. Deviations of the observed diversity from the 
null expectation indicated a nonrandom spatial distribution of 
fungi or bacteria at a given scale. 

Results 

Soil property variations 

Soil pH comparisons using a Mann-Whitney U test of 
significance, between five plots of natural and five plots of 
converted TRF ecosystems in each of the four forest sites, found 
significantiy higher (less acidic) pH in three of the four sites at 
Budongo (/) = 0.0107), Kaniyo Pabidi (/> = 0.0112), and Mabira 
ip = 0.0269); and non-significant difference at Maramagambo 
(jf) = 0.1706). Percentage soil organic carbon (SOC) was signifi- 

Inventorv diversity Differential diversity Relevant scale 

Qi (vyithin plots) + t Pi (Between plots) plots (samples) 

1- 

a2 (within sites) P2 (between sites) Sites 

1- 

aj (within ecosystem type) P3 (between ecosystem types) Ecosystem types 
Y (total diversity at ecosystem type level)^ l 

Figure 2. Illustration of hierarchical spatial scales in our 
additive partitioning model. The a scale is the within-level, and |3 
scale, the between-level components. Because a diversity at a given 
scale is the sum of the a and (3 diversity at the next lower scale, the total 
diversity (y) can be described by the following formula: 0(1+pi+p2+p3 
[14,22]. 

doi:1 0.1 371 /journal.pone.Ol 0481 8.g002 



cantly higher in natural than converted ecosystems in all four sites 
i.e. Budongo (/) = 0.0119), Kaniyo Pabidi (/) = 0.0212), Mabira 
(/j = 0.0122) and Maramagambo lp = 0.0119) with combined 
%SOC in natural sites more than double of that in converted 
sites; whereas %soil nitrogen was only significantly higher in 
natural forests at Budongo (^ = 0.0112) and Kaniyo Pabidi 
(/) = 0.0119), and non-significant at Mabira (^ = 0.6015) and 
Maramagambo (p = 0.0947) (Table 3). 

Ecosystem and site connparisons of microbial community 
diversity 

Bacterial (B) communities were significantly richer (p = 0.0304; 
Mann- Whitney U) in detectable bands than fungal (F) communi- 
ties in both converted (C) and natural (N) ecosystems (converted: 
medians; F = 36, B = 61.5; natural: medians; F = 39.5, B = 60.5). 
While total band richness (B-hF) did not differ between natural and 
converted forests we observed greater fungal richness in natural 
than in converted forests (medians: C = 36, N = 39.5; test 
Stat = 18.5) and more bacterial bands in converted than in natural 
ecosystems (medians: C = 61.5, N = 60.5; test stat= 18.5). Kaniyo 
Pabidi was the most diverse site overall with the highest number of 
bacterial and fungal bands, while Maramagambo had the least 
band richness (Fig 3). 

Natural sites harbored more bands unique to one site than 
converted sites for bacteria at Kaniyo Pabidi and Maramagambo 
and for fungi at Maramagambo and Budongo. Mabira and 
Kaniyo Pabidi had higher numbers of unique bacterial bands than 
at Maramagambo and Budongo. There were also more unique 
fungal bands at Mabira and Budongo than at Maramagambo and 
Kaniyo Pabidi (Fig 4). 

We also found that Mabira and Maramagambo had the lowest 
bacterial Jaccard's community similarity indices [32] between 
natural and converted ecosystems, whereas Budongo and Mabira 
had the lowest fungal community similarity between natural and 
converted ecosystems (Table 4). Dissimilarity between natural and 
converted ecosystems was nonetheless non-significant in all sites 
for both fungal and bacterial communities. Also, there was 
generally greater dissimilarity between sites of fungal communities 
than in bacterial communities suggesting a higher susceptibility to 
habitat change among fungi than bacteria (Table 4). 

Ecosystem classification and importance of predictor 
variables 

The CCA showed that, despite the relatively small amount of 
difference between sites, soil pH, average phosphorus, and texture 
(%sand) had strong influence on bacterial diversity in the TRF 
ecosystem (Fig 5); whereas organic carbon, sodium, pH and 
average phosphorus were strongly associated with fungal commu- 
nity variation in both natural and converted TRF ecosystems 
(Fig 5). The CCA also showed that bacterial communities in both 
Kaniyo Pabidi and Mabira were unique to bacterial communities 
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in the other sites and there was high contrast between bacterial 
communities of converted and natural ecosystems at Kaniyo 
Pabidi. Fungal communities at Maramagambo and Mabira were 
also unique to those in other sites and there was high contrast 
between fungal communities at Mabira's natural and converted 
ecosystems. Furthermore, the CCA showed that fungal commu- 
nities responded more to soil pH levels than bacterial communities 
(Fig 5), with site-specific patterns showing that bacteria and fungi 
were grouping according to sites. 

A discriminant analysis to predict whether bacterial or fungal 
communities were from natural or converted ecosystems found 
that only OC, Ca, N, and pH for bacterial communities; and OC, 
N, Ca, and pH for fungal communities (all ranked from most 
important to least important) were found to be significant 
predictors of soil physicochemical properties. AH other variables 
were poor predictors in this context (Table 5). 

Hierarchical scaling 

We found 58 and 56 fungal bands in natural and converted 
forests respectively, from 1 7 plots of natural ecosystems and 20 
plots of converted forests. There were also 92 and 88 bacterial 
bands in natural and converted ecosystems respectively found in 
20 plots of converted ecosystems and 17 plots of natural ecosystem. 
All these were within four sites, ^-diversity varied more than ot- 
diversity between natural and converted ecosystems for both 
bacteria and fungi. We found higher bacterial and fungal (3- 
diversity in converted ecosystems than in natural ecosystems at 
lower hierarchical scales (Pi); higher P-diversity in natural than 
converted at between-site scale (P2), and higher P-diversity in 
converted than in natural ecosystems at the between-ecosystem 
type scale (P3) (Fig 6). 

We also found substantial contribution of observed P-diversity 
(Pi, P2, and Pa) to total band richness (y-diversity), while a- 
diversity of both bacteria and fungi in converted and natural 
ecosystems were similar. Spatial partitioning of total diversity also 
consistently showed that the beta components (Pi and P2) were 
always greater than expected by chance, whereas the alpha 
component (oti) was always lower than expected. For both fungal 
and bacterial communities in natural and converted ecosystems, 
observed within plot diversity were substantially less than values 
expected from individual-based randomizations (Fig 7). 

Discussion 

Soil property variations and site differences 

Studies in both tropical and temperate zones show that soils 
in converted or cropped areas normally have reduced soil 
aggregation, structural stability and organic matter, and an 
increase in bulk density when compared to forests [40,41]. 
Habitat conversion may also alter soil properties such as 
nutrient levels, and abiotic conditions and may affect associ- 
ations between organisms. In our study there are some local 
details that may influence our results. 

Both Maramagambo and Kaniyo Pabidi are located within 
Queen Elizabeth NP and Murchison Falls NP respectively and 
are protected by Uganda Wildlife Authority (UWA) personnel. 
They are well protected and there is little evidence of recent 
encroachment. There is significant wildlife populations includ- 
ing elephants, buffaloes, zebras and the areas are frequented 
by tourists. Protection by UWA and presence of dangerous 
animals (such as buffalos and lions) reduce damaging human 
activity at Kaniyo Pabidi and Maramagambo which should 
enhance the difference between natural and converted 
ecosystems. Maramagambo's location, in contrast, means the 
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■ Fungi, C 
Q Fungi, N 

■ Bacteria, C 
□ Bacteria, N 



Mabira Maramagambo Budongo Kaniyo Pabidi 



Figure 3. Band richness for fungal and bacterial communities in converted and natural ecosystems. All richness values are total bands 

present in five samples of each ecosystem treatment (error bars are 5% confidence interval). 
doi:1 0.1 371/journal.pone.01 0481 8.g003 



forest is subjected to substantial natural disturbance from 
frequent storms and strong erosive runoff even within the 
natural forest, whereas tourist activity at Kaniyo Pabidi 
seemed to have little impact on soil properties. Converted 
areas at Kaniyo Pabidi were also sparsely populated with 
limited human impacts on the environment. Its sites were old 
and might have been cultivated for at least 20years. 



In our study, Budongo is located next to a high, mainly 
subsistence population and resultant population activity. But even 
though encroachment, illegal hunting and logging in natural 
habitats in Budongo are not uncommon, there seems to be 
minimal impact of conversion on soil properties in our sample 
locations; whereas proximity of Mabira's natural forest to densely 
populated urban areas exposes it to increased human activities, 
likely reducing its difference with converted sites. 




Figure 4. Bacterial and fungal bands unique to converted (C) and natural (N) ecosystems at each site (error bars are 5% confidence 
interval). 

doi:1 0.1 371/journal.pone.010481 8.g004 
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Microbial community variations 

Soil properties determine many aspects of soil microbial 
community structure [42-44]. Carbon availability [45-47], 
nitrogen availability [45,48,49] and soil pH [44,50,51] can all 
influence microbial community composition and diversity. In 
addition, correlation studies have shown that plant species [52-54] 
and soil type [45,54,55] are associated with variation in microbial 
communities. It has also been shown that land use indirectiy affects 
bacterial community structure by modification of soil properties 
[56] but similarity between converted and natural ecosystem 
bacterial communities may also suggest a high number of 
generalists. 

Nacke et al., [56] found that bacterial commvmity composition 
in forests and grassland was largely determined by tree species and 
soil pH. Jesus et al. [6] also showed that bacterial community 
structure is influenced by changes linked to soil acidity and 
nutrient concentration. Other studies also suggest that soil pH is a 
major factor influencing microbial community composition 
[50,57—59]. This influence of soil pH has been recognized at 
different taxonomic levels [45,60] with most microorganisms 
thriving within a limited pH range. This is because acids can 
denature proteins and large pH changes may inhibit growth in 
microorganisms. Fierer and Jackson, 2006 [44] found, in contrast, 
that net carbon mineralization rate (an index of C availability) was 
the best predictor of phylum-level abundances of dominant 
bacterial groups, and Bisset et al., [42] found that soil microbial 
communities were consistent with disturbance gradients within 
different agricultural treatments and relatively undisturbed non- 
agricultural sites. 

Because of widespread forest conversion in Uganda as a 
result of increasing population pressure, estimated at between 
1.1% and 3.15% per year [61]; natural ecosystems and their 
inhabitant biodiversity are at risk [62]. Loss of diversity 
increases the likelihood of losing important functional roles 
and associated ecosystem processes. At landscape scale, spatial 
and temporal variations of microbial communities in forest 
soils are influenced by numerous biotic and abiotic factors. 
These factors may include climate, soil types, and vegetation 
associations [50,63,64]. Owing to this study design, many of 
these factors were assumed to be similar between natural and 
converted ecosystems. For instance, the proximity of natural 
and converted ecosystem sites meant that climate and geology 
were, we assume, similar in the two treatments. Even though 
there could still be a number of underlying causes of 
community differences, two likely influences were assumed to 
be soU properties [45,54,55] and vegetation types [52-54]. 
Clearly both of these sets of factors change when forest is 
converted for agriculture or range lands. 

Despite substantial reductions in SOC, N, Ca and pH in 
converted sites in this study, differences in microbial communities 
were small meaning that converted sites still had sufficient SOC, 
N and Ca to sustain the same microbial populations. The close 
proximity of the matched pairs could also lead to a source-sink 
relationship between the natural and converted forests, with the 
presence of unique taxonomic groups a likely indication of 
habitat preference (endemism) for some taxonomic groups. It 
may also be an indicator of relative habitat dissimilarity. The high 
numbers of unique bacterial bands at Mabira and Kaniyo Pabidi 
and unique fungal bands at Mabira and Budongo (Fig 4) thus 
suggests that ecosystem alteration at these sites was sufficient to 
force a different regime of processes and structures enabling a 
new set of taxonomic groups to predominate. Mabira had high 
numbers of both unique bacterial and unique fungal bands that 
can be attributed to the extent of disturbance at its sites (Mabira 
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Figure 5. CCA for bacterial (B) and fungal (F) relationships using relative abundance of bands and soil physicochemical properties 
in natural and converted ecosystems. The symbols (left graphs) represent the similarity between each sample (plot) as defined by their diversity, 
and the vectors (right graphs) represent the structural matrix for soil properties and their influence on relative abundance of each band. The length of 
the vectors represents the relative strength of influence of the particular aspect of soil physicochemical property. 
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is the only peri-urban tropical rainforest site among the four 
selected sites). 

The low numbers of unique bacterial and fungal bands at 
Maramagambo can be attributed to the high erosion at natural 
sites that reduced the contrast between the natural and converted 
sites. For the other sites, bacteria and fungi had different responses 
to ecosystem alterations. This could indicate separate influences on 
microbial distribution that exist when alteration is moderate. 
Similarity indices suggested that bacterial and fungal communities 
were determined by separate forces leading to cUstinct responses 
across the study locations. 

Hierarchical scaling 

Many studies have shown that specialist species are more 
negatively affected by current global changes than generalists 
[7,65]. The process of biotic homogenization can involve the 
replacement of native biota with non-natives or the introduction of 
generalist species [66] . In this study, the net decrease in fi-diversity 
from natural to converted TRF ecosystem at the between-site scale 
(fi2) for both fungi and bacteria was an indication of biotic 
homogenization [18,66]. This can result from ecosystem alter- 
ations which can in-turn alter ecosystem function and reduce 
ecosystem resilience to disturbance [65,67]. 



We also showed that the fi components of diversity (Pi and P2; 
the average diversity between the plots and sites, respectively) were 
consistently higher than those expected by chance, whereas the 
local «! diversity component (ofi, the average diversity within the 
plots) was consistently lower than that expected (Fig 7). Such scale- 
dependent deviations of the observed diversity from the expected 
can be generally explained by aggregation at a relatively small 
"local" scale and, spatial dififerentiation of diversity at a larger 
"landscape" scale [33,34,68,69]. 

Relatively lower diversity within converted ecosystems suggests 
that conversion of natural TRF ecosystems results in reduced 
diversity for both bacteria and fungi. This is compatible with 
recent studies that show that conversion of TRF ecosystems 
threatens microbial diversity [7] and because microorganisms, like 
all other organisms, have habitat preferences and may be affected 
by land-use changes [6,64]. While we cannot be certain that such 
decline in diversity has led to a decline in any particular ecosystem 
functions or services, this is a possibility that deserves further 
evaluation, and we speculate that such loss of diversity will at the 
very least cause a reduction in functional redundancy and 
associated resilience. 

Higher P-diversity of both bacterial and fungal communities at 
the between-plot scale (Pi) in converted ecosystems indicates 
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Table 5. Structure matrix rank showing absolute size of correlation between discriminant analysis function from most important 
to least important predictor variable of site factors (soil physiochemical properties) and their influence on the variation of soil 
microbial communities. 



Bacteria 




Fungi 




Predictor Variables 


Function 1 


Predictor variables 


Function 1 


OC 


0.625* 


OC 


0.625* 


Ca 


0.471* 


N 


0.493* 
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0.421* 


Ca 


0.473* 


pH 


0.355* 


pH 


0.340* 


Mg 


0.298 


Mg 


0.281 


Cu 


0.197 


Cu 


0.221 


Na 


0.181 


Na 


0.183 
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0.169 


K 


0.160 


Av.P 


-0.077 


Sand 


0.064 


Sand 


0.070 


Av.P 


-0.060 


Simpson 


-0.065 


Fe 


0.024 


Fe 


0.029 


Shannon 


0.012 


Shannon 


-0.023 


Simpson 


0.011 



(* = important predictor variable, with 0.30 used as the threshold). 
doi:l 0.1 371 /journal.pone.Ol 0481 8.t005 



difTerentiation (reduced community similarity) in converted determining total richness and community composition [14,34]. 
ecosystems at this hierarchical scale. Considering the multiple Differences between the observed and expected diversity compo- 
land-uses and cropping systems of converted areas, this was nents could be due to ecological processes that lead to a non- 
expected. There was also substantial contribution of ^-diversity to random dispersion of individuals. These processes could include 
total diversity (y). This suggests the importance of nonrandom intra-specific aggregation, habitat selection, and limited dispersal 
ecological processes at the between-plot and between-site scale in capacity [33] . 



Alpha diversity for bacterial Beta diversity for bacterial communities 




Sampling levels Sampling levels 

Figure 6. Additive partitioning of bacterial and fungal diversity (expressed as additive richness) across alpha and beta hierarchical 
spatial scales at three sampling levels (plot, site and ecosystem type) in natural and converted TRF ecosystems. 

doi:1 0.1 371/journal.pone.Ol 0481 8.g006 
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Figure 7. The additive partitioning of total bacterial and fungal community, y -diversity into a and p components at three nested 
spatial scales, with each component expressing their relative contributions to y -diversity; where y -diversity is equal to 

a1+pi+p2+p3. The observed (obs) partitions are compared with the expected (exp) values, as predicted by the null model based on 1000 iterations 
using individual-based randomization. 
doi:l 0.1 371/journal.pone.01 0481 8.g007 



Conclusion 

There is international concern about the threat to natural 
habitats in the Nile river watershed and the consequential loss of 
important biodiversity. Whereas aspects of microbial biogeogra- 
phy and influence of forest conversion in Uganda's Nile river 
watershed is largely unknown, this study offers an important first 
glimpse into indicators of spatial diversity patterns of soil fungal 
and bacterial communities in the Uganda's Nile river watershed. 
Our observations of reduced soil microbial diversity, both bacterial 
and fungal, in converted ecosystems though unsurprising in itself 
causes us some concern and would justify further work to 
determine the significance of the diversity lost and the wider 
implications. 

By focusing on diversity patterns across multiple hierarchical 
spatial scales, we were able to identify the scale at which 
regional microbial diversity is maximized. We showed that there 
was substantial contribution of P-diversity to total ecosystem 
diversity (y) which includes taxa at the between-plot, site and 
ecosystem scales and unique taxa, highlighting the necessity to 
conserve marginal habitats and ecotones. Soil microbial 
communities in Uganda's Nile river watershed exhibit consid- 
erable resilience to forest conversion even though SOC, N, Ca 
and pH were all significantly altered. This result is surprising 
given that these physical and chemical properties typically 
strongly influence microbial diversity. Additionally, the varia- 
tion among sites was quite large, indicating that soil commu- 
nities in this region vary considerably on a regional spatial scale. 
Our results do not explain this variation. Most studies suggest 
that biogeographic barriers play little role in the geographic 
structure of soil communities. Rather than a consistent general 
pattern of microbial community change following forest 
conversion we find that responses are largely site-specific and 
widely variable. 
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